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Abstract 

We investigate the finite temperature expectation values of tlie cliarge and current densities 
for a complex scalar field with nonzero chemical potential in background of a flat spacetime with 
spatial topology RF x {S^Y. Along compact dimensions quasiperiodicity conditions with general 
phases are imposed on the field. In addition, we assume the presence of a constant gauge field 
which, due to the nontrivial topology of background space, leads to Aharonov-Bohm-like effects 
on the expectation values. By using the Abel-Plana-type summation formula and zeta function 
techniques, two different representations are provided for both the current and charge densities. 
The current density has nonzero components along the compact dimensions only and, in the absence 
of a gauge field, it vanishes for special cases of twisted and untwisted scalar fields. In the high- 
temperature limit, the current density and the topological part in the charge density are linear 
functions of the temperature. The Bose-Einstein condensation for a fixed value of the charge is 
discussed. The expression for the chemical potential is given in terms of the lengths of compact 
dimensions, temperature and gauge field. It is shown that the parameters of the phase transition 
can be controlled by tuning the gauge field. The separate contributions to the charge and current 
densities coming from the Bose-Einstein condensate and from excited states are also investigated. 

PACS numbers: 03.70. -hk, ll.lO.Kk, 03.75.Hh 



1 Introduction 
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In recent years, there has been a large interest to the physical problems with compact spatial dimen- 
sions. Several models of this sort appear in high energy physics, in cosmology and in condensed matter 
physics. In particular, many of high energy theories of fundamental physics, including supergravity 
and superstring theories, are formulated in spacetimes having extra compact dimensions which are 
characterized by extremely small length scales. These theories provide an attractive framework for the 
unification of gravitational and gauge interactions. The models of a compact universe with nontrivial 
topology may also play an important role by providing proper initial conditions for inflation p!]. 

In the models with compact dimensions, the nontrivial topology of background space can have 
important physical implications in classical and quantum field theories, which include instabilities in 
interacting field theories [2j , topological mass generation [U |4] and symmetry breaking [U [5] . The 
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periodicity conditions imposed on fields along compact dimensions allow only the normal modes with 
suitable wavelengths. As a result of this, the expectation values of various physical observables are 
modified. In particular, many authors have investigated the effects of vacuum or Casimir energies 
and stresses associated with the presence of compact dimensions (for reviews see Refs. [6], [?])• The 
topological Casimir effect is a physical example of the connection between quantum phenomena and 
global properties of spacetime. The Casimir energy of bulk fields induces a non-trivial potential for the 
compactification radius of higher-dimensional field theories providing a stabilization mechanism for 
the corresponding moduli fields and thereby fixing the effective gauge couplings. The Casimir effect 
has also been considered as a possible origin for the dark energy in both Kaluza-Klein type models 
and in braneworld scenario [8]. 

The main part of the papers, devoted to the infiuence of the nontrivial topology on the properties 
of the quantum vacuum, considers the vacuum energy and stresses. These quantities are chosen 
because of their close connection with the structure of spacetime through the theory of gravitation. 
For charged fields another important characteristic, bilinear in the field, is the expectation value of 
the current density in a given state. In Ref. [9], we have investigated the vacuum expectation value of 
the current density for a fermionic field in spaces with an arbitrary number of toroidally compactified 
dimensions. Application of the general results are given to the electrons of a graphene sheet rolled into 
cylindrical and toroidal shapes. For the description of the relevant low-energy degrees of freedom we 
have used the effective field theory treatment of graphene in terms of a pair of Dirac fermions. For this 
model one has the topologies x and (S*^)^ for cylindrical and toroidal nanotubes respectively. 
Combined effects of compact spatial dimensions and boundaries on the vacuum expectation values 
of the fermionic current have been discussed recently in Ref. [lOj . In the latter, the geometry of 
boundaries is given by two parallel plates on which the fermion field obeys bag boundary conditions. 
The effects of nontrivial topology around a conical defect on the current induced by a magnetic flux 
were investigated in Ref. pULj for scalar and fermion flelds. 

In the present paper we consider the finite temperature charge and current densities for a scalar 
field in background spacetime with spatial topology RP x {S^y. In both types of models with compact 
dimensions used in the cosmology of the early Universe and in condensed matter physics, the effects 
induced by the finite temperature play an important role. The thermal corrections arise from thermal 
excitations of fluctuation spectrum and they depend strongly on the geometry. As a consequence of 
this, thermal modifications of quantum topological effects can differ qualitatively for different geome- 
tries. The thermal Casimir effect in cosmological models with nontrivial topology has been considered 
in Refs. [12]. A general discussion of the finite temperature effects for a scalar field in higher di- 
mensional product manifolds with compact subspaces is given in Ref. [13]. Specific calculations 
are presented for the cases when the internal space is a torus or a sphere. In Ref. [14j . the corre- 
sponding results are extended to the case in which a chemical potential is present. In the previous 
discussions of the effects from nontrivial topology and finite temperature, the authors mainly consider 
periodicity and antiperiodicity conditions imposed on the field along compact dimensions. The latter 
correspond to untwisted and twisted configurations of fields respectively. In this case the current 
density corresponding to a conserved charge associated with an internal symmetry vanishes. As it will 
be seen below, the presence of a constant gauge field, interacting with a charged quantum field, will 
induce a nontrivial phase in the periodicity conditions along compact dimensions. As a consequence 
of this, nonzero components of the current density appear along compact dimensions. This is a sort 
of Aharonov-Bohm-like effect related to the nontrivial topology of the background space. 

The organization of the paper is as follows. In the next section the geometry of the problem 
is described and the thermal Hadamard function is evaluated for a complex scalar field in thermal 
equilibrium. In Section [3l by using the expression for the Hadamard function, the expectation values 
of the charge and current densities are investigated. Various limiting cases are discussed. Alternative 
expressions for the charge and current densities are provided in Section [5] by making use of the zeta 
function renormalization approach. The Section [6] is devoted to the investigation of the Bose-Einstein 
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condensation in the background under consideration. The properties of the vacuum expectation 
value of the charge density are discussed in Appendix El Throughout the paper we use the units 
h = c = kB = 1, with ks been Boltzmann constant. 

2 Geometry of the problem and the Hadamard function 

We consider the quantum scalar field (p{x) on background of {D + 1) dimensional flat spacetime 
with spatial topology i?^ x (S^)'^ ,p + q = D. For the Cartesian coordinates along uncompactified and 
compactified dimensions we use the notations Xp = {x^,...,xP) and Xg = (x^^^, x^), respectively. 
The length of the l—ih compact dimension we denote as Li. Hence, for coordinates one has — oo < 
< c« for / = 1, and Q ^ ^ Li for / = p + 1, ...,D. In the presence of a gauge field the 
field equation has the form 

{g^^'D^D, + m^)ip = 0, (2.1) 

where = 5^ + ieA^^ and e is the charge associated with the field. One of the characteristic features 
of field theory on backgrounds with nontrivial topology is the appearance of topologically inequivalent 
field configurations [15j. The boundary conditions should be specified along the compact dimensions 
for the theory to be defined. We assume that the field obeys generic quasiperiodic boundary conditions, 

(/7(t,Xp,Xg + Liei) = e*"V(*>Xp,Xg), (2.2) 

with constant phases \ai\ ^ vr and with being the unit vector along the direction of the coordinate 
x\l = p + 1, ...,D. The condition (|2.2p includes the periodicity conditions for both untwisted and 
twisted scalar fields as special cases with a; = and a; = vr, respectively. 

In the discussion below we will assume a constant gauge field A^. Though the corresponding 
field strength vanishes, the nontrivial topology of the background spacetime leads to the Aharonov- 
Bohm-like effects on physical observables. In the case of constant A^, by making use of the gauge 
transformation 

^{x) = e--V'(^), A^ = A'^ + d^x, (2.3) 

with X = Afj^x^ we see that in the new gauge one has A'^ = and the vector potential disappears from 
the equation for ip'{x). For the new field we have the periodicity condition 

(f'it,xp,xq + Liei) = e'^'^'{t,xp,xg), (2.4) 

where 

ai = ai + eAiLi. (2.5) 

In what follows we will work with the field <^'{x) omitting the prime. Note that for this field = 3^. 
As it is seen from Eq. ()2.5p . the presence of a constant gauge field shifts the phases in the periodicity 
conditions along compact dimensions. In particular, a nontrivial phase is induced for special cases 
of twisted and untwisted fields. As it will be shown below, this is crucial for the appearance of the 
nonzero current density along compact dimensions. Note that the term in Eq. (j2.5p due to the gauge 
field may be written as 

eAiLi = 2tt^i/^o, (2.6) 

where is a formal flux enclosed by the circle corresponding to the l-th. compact dimension and 
$0 = 2vr/e is the flux quantum. 

The complete set of positive- and negative-energy solutions for the problem under consideration 
can be written in the form of plane waves: 

^t\x) = C^e'''-'^'^\ iOk = + m2, (2.7) 
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where k = (kp, kg), kp = (fci, . . . , kp), kg = (/cp+i, . . . , A:^)), with — oo < ki < +00 for i = 1, . . . ,p. For 
the momentum components along the compact dimensions the eigenvalues are determined from the 
conditions (12. 4|) : 

ki = {2TTni + di)/Li, ni = 0,±l,±2,...., (2.8) 



with I = p + 1, D. From Eq. (j2.8p it follows that the physical results will depend on the fractional 
part of Qi/(27r) only. The integer part can be absorbed by the redefinition of rii. Hence, without loss 
of generality, we can assume that \ai\ ^ vr. The normalization coefficient in ()2.7p is found from the 
orthonormalization condition 

d''x^l'\x)^^^'>ix) = ^-^AA'^k', (2.9) 

where J^k' = S{kp — kp)5„^^^ Substituting the functions (j2.7p . for the normalization 

coefficient we find 

with Vq = Lp^i....L£) being the volume of the compact subspace and 

= yi^^T^^, k^ = x; ('-^^y. (2.11) 

l=p+l ^ '■ ^ 

The smallest value for the energy we will denote by wq- Assuming that la^l ^ vr, we have 



^o = A/Ef=p+i«'/^'+"^'- (2.12) 

We are interested in the expectation values of the charge and current densities for the field (fix) in 
thermal equilibrium at finite temperature T . These quantities can be evaluated by using the thermal 
Hadamard function 

GW(x,x') = ((/7(x)vp+(x') + (/'+(x')</'(x)> 

= tr[;5((^(x)(p+(x')+V'+(x')</^(x))], (2.13) 

where (• • • ) means the ensemble average and p is the density matrix. For the thermodynamical 
equilibrium distribution at temperature T, the latter is given by 

p = Z-ie-^(^-'^'«), (2.14) 

where /3 = 1/T. In Eq. ()2.14p . Q denotes a conserved charge, /i' is the related chemical potential and 
Z is the grand-canonical partition function 

Z = tr[e-^(^-^'<3)]. (2.15) 

In order to evaluate the expectation value in Eq. (|2.13p we expand the field operator over a 
complete set of solutions: 

ip{x) = ^[a^yp'^\x) + b+ipi-\x)], (2.16) 

k 

with Ylk = / d\<.p X^n, ~ ("-p+i' • • • 1 ^d)- Here and in what follows we use the notation 

+00 +00 

n 711 = — 00 71; = — 00 
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for n = (ni, . . . ,ni). Substituting the expansion (|2.16p into Eq. (|2.13|) . we use the relations 

where = efi' . Note that the chemical potentials have opposite signs for particles (^u) and antiparticles 
(— /x). The expectation values for Skaj^; and bi^b^, are obtained from (j2.18p by using the commutation 
relations and the expectation values for the other products are zero. For the Hadamard function we 
get 

0« = OW + E i^^S^ , (2.19) 

k s=± 

where the first term in the right-hand side corresponds to the zero temperature Hadamard function: 



G^f^\x,x') = {0\ip{x)ip+{x') + ip+{x')ip{x)\0) 

= EE^k^(^)^^^*(^'), (2.20) 

k s=± 



with |0) being the vacuum state. In order to ensure a positive-definite value for the number of particles 
we assume that ^ wq, where ujq is the smallest value of the energy (see Eq. (j2.12p ). 

By using the expressions (|2.7p for the mode functions and the expansion (e^ — = Yl'n=i 
the mode sum for the Hadamard function is written in the form 

UH-n ih- .A^ \ ^ c ^ 



^ ' ^ Vj {27t)p ^ 



cos(LJkAt) + X] ^' 



1 



^uJi^{siAt—nl3)—snfif} 

n=l s=± 



(2.21) 



where Axp= Xp — x^, Axg= Xg — x^, At = t — t' . For the evaluation of the Hadamard function we 
apply to the series over rir the Abel-Plana-type summation formula |16t [T7] (for applications of the 
Abel-Plana formula and its generalizations in quantum field theory see Refs. [6l 118^ [T9]) 



2 °° fOO 

— J2 9{kr)f{\kr\)= dz[g{z)+g{-z)]f{2 
^ ^ — — , 



nr=— oo 

f CXD 



+^j dz[/(zz)-/(-z.)] 5: (2.22) 



pZLr+iXctr 1 ' 

A=±l 



where kr is given by Eq. (j2.8p . For the Hadamard function we find the expression 



G'"(-^') = G<V,(...') + ^/(|^E 



gikp-Axp+ikq_i-Ax,_i 

9-1 



n=— oo 



~ .oo cosh[(At-m/3)Jz2-a;2 ] -\zAx- 
X y e'^"'^ / dz ^ — - y -, (2.23) 



where n^_^ = (n^+i, . . . , n^_i, n,.+i, . . . , n^), k<^-i = (fcp+i, . . . , fc^-i, A:r+i, • • • , /^d), and 



-p,,_l = Jk2+k2 +^2. (2.24) 
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The first term in tlie right-liand side of Eq. (j2.23p . G^^^ q_^{x,x'), comes from the first term on the 
right of Eq. (j2.22p and it is the Hadamard function for the topology RP^^ x (S"^)^^^ with the lengths 
of the compact dimensions (ip+i, . . . , Lr-i, Lr+i, . . . , Ld). 

For the further transformation of the expression ()2.23p we use the expansion 



_ _ ^-z{lLr+XAx''')-\ilar ^2 25) 



1=1 

With this expansion the z-integral is expressed in terms of the Macdonald function of the zeroth order. 
Then the integral over kp is evaluated by using the formula from Ref . [20j . For the Hadamard function 
we arrive to the final expression 

2L °° 

C^^^(x x') = ^ ^ \^ gm,.Or+ra/i/3gikg_i-Ax,_i 

' (27r)P/2+l ^ ^ 

^ ' n=— 00 rig 

X^lr_Jp/2 ^\^^P? + (Ax'^ - UrLrf - {At - in(3f^ , (2.26) 

where 



Mx) = x-''K,{x), ujnr_^ = ^k2_, + m2. (2.27) 

Note that the = term in Eq. (j2.26p corresponds to the function G^^^ ^_^{x,x'). Hence, the part 
of the Hadamard function in Eq. (|2.26p with n,. 7^ is induced by the compactification of the r-th 
direction to a circle with the length L^. 

An alternative expression for the Hadamard function is obtained directly from Eq. (|2.2ip . We 
first integrate over the angular part of kp and then the integral over |kp| is expressed in terms of the 
Macdonald function. The corresponding expression is written in terms of the function (I2.27P as 



(2^)- 



x/p^(u;n,A/|A^pl' - - (2-28) 



uJn, = \K + m2, (2.29) 



with the notation 



and kg is given by Eq. (j2.1ip . Note that the explicit information contained in Eq. (|2.26p is more 
detailed. Both representations (12.26P and (I2.28P present the thermal Hadamard function as an infinite 
imaginary-time image sum of the zero temperature Hadamard function. This is the well-known result 
in finite temperature field theory (see, for instance, Ref. |21j). 

3 Charge density 

Having the thermal Hadamard function we can evaluate the expectation value for the current density 

ji{x) = ie[ip^ {x)diip{x) - (5/93+(x))99(x)], (3.1) 
/ = 0, 1, . . . , D, by using the formula 

{ji{x))='-e]im{di-d[)G^^\x,x'). (3.2) 
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By making use of the relation dzfy{z) = —zf^^i{z), from Eq. (|2.26|) for the charge density (/ = 0) 
one finds 

Oo) = -— -p-f— V cos(n^Q^) Vnsinh(n/i/3) 

X E ^n^Vf +1K-1 Vn'rL'r + n^n, (3.3) 

where the prime on the sign of sum means that the term = should be taken with the coefficient 
1/2. 

As it is seen from Eq. (j3.3|) . the charge density is an even function of the phases q; and, for a fixed 
value of the chemical potential, it vanishes in the zero temperature limit. It is a periodic function of 
di with the period equal to 2ir. In the case of zero chemical potential the charge density is zero. In 
Eq. (jS.Sp . the term with n,, = corresponds to the charge density for the topology R^^^ x (S^)'^^^ 
with the lengths of the compact dimensions (i^p+i, • • • , Lr-i, Lr+i, . . . , L/j) and the contribution of 
the terms with rir 7^ is the change in the charge density due to the compactification of the r-th 
dimension to with the length L^. By taking into account Eq. (12.6(1 . we see that the charge density 
is a periodic function of fluxes with the period equal to the flux quantum. Note that the sign of 
the ratio (j'o) / e coincides with the sign of the chemical potential. 

An alternative expression for the charge density, more symmetric with respect to the compact 
dimensions, is obtained by applying the formula 

+00 ^ +00 



J2 cos{na)McV^Ta^) = ^ E ^n~'fu~i/2{bwn), (3.4) 



n=— 00 



with a, 6, c > 0, Wn = \/ {2-Kn + aY ja^ + c^, to the series over in Eq. ()3.3p . This leads to the 
expression 

(io) = 5]nsinh(n/./3) ^<+Vrii(n/3u;„,), (3.5) 

(27r) 2 „=i 

with the notation (j2.29p . This formula could also be directly obtained from Eq. (j3.2p using the 
expression (|2.28p for the Hadamard function. The form (|3.5p for the charge density in the case of 
topology X {S^yi~^ is also obtained from Eq. ()3.3p taking the limit — )• 00. 

In the case of Minkowski spacetime one has p = 9 = 0, and from Eq. (j3.5p we get 

Oo)(M) = V?^sinh(n/i/3)/D+i(n/3m), (3.6) 

(27r)— 

with \^\ ^ m. The thermodynamic properties of the relativistic Bose gas in this case have been 
considered in Refs. [22J 123] . If all spatial dimensions are compactified, the corresponding formulas are 
obtained from Eqs. (13. 3p and (13. 5p taking p = 0. In particular, from Eq. (|3.5p one has 

00 

Oo) = T7 E sinh(n^/3) ^ e-^P^-^ , (3.7) 

^1 n=l 

where we have used fi/2{x) = y7r/2x~^e~^. 

Let us consider some limiting cases of Eq. ()3.5p . If the length of the /-th compact dimensions is 
large compared to other length scales, in the sum over ni in Eq. (j3.5p the contribution from large 
values of ni dominates and, to the leading order, we replace the summation by the integration. The 
corresponding integral is evaluated with the help of the formula 

dy{y'' + b''Y-^ f,_^{c^^J^) = ^lf+^f._^M), (3.8) 
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and from Eq. (|3.5|) we obtain the expression of the charge density for the topology R^^^ x (S^y^^. 

If the length of the l-th. compact dimension is small compared with the other length scales and 
Li <^ P, under the assumption \ai\ < vr, the main contribution to the corresponding series in Eq. (13. Sh 
comes from the term with ni = 0. The behavior of the charge density is essentially different with 
dependence whether the phase ai is zero or not. When ai = 0, we can see that, to the leading order, 
Li {jo) coincides with the charge density in [D — l)-dimensional space of topology x (S^)'^^^ and 
with the lengths of the compact dimensions Lp^i,. . . ,Li_i,L/+i,. . . ,L£). In particular, this is the case 
for an untwisted scalar field in the absence of a gauge field. For d/ 7^ and for small values of L/, the 
argument of the Macdonald function in Eq. (13. Sp is large and the charge density is suppressed by the 
factor e^l"'!'^/^'. 

In the low-temperature limit the parameter /3 is large and the dominant contribution to the charge 
density comes from the term n = 1 in the series over n and from the term in the series over with 
the smallest value of which corresponds to n/ = 0, / = p + 1, . . . , £). To the leading order we find 

(2vr)p/2+i/3p/2'^o ' ^ ^ 

with loq given by Eq. (j2.12p . 

From Eq. (j3.5p it follows that the expectation value of the charge density is finite in the limit 
— )• ujo for p > 2 and it diverges for p ^ 2. In order to find the asymptotic behavior near the 
point |;u| = loq, we note that for p ^ 2, under the condition /3{uJo — \p\) <^ 1, the main contribution 
to Eq. (|3.5p comes from the term with n; = (wn, = ^0) and in the corresponding series over n 
the contribution of large n dominates. In this case we can use the asymptotic expression for the 
Macdonald function for large values of the argument and to the leading order this gives: 

Oo> « sgn(M):^ [^j Li,/2(e-/^(-o-l'^l)), (3.10) 



where Lis(x) is the polylogarithm function. For the latter one has Lio(x) = x/{l — x), Lii(2;) = 
— ln(l — x). By taking into account that Lis(e~^) w r(l — s)y'^~^ for \y\ <C 1 and s < 1, one finds the 
following asymptotic expressions: 

^ sgn(/i)r(l-p/2) fu;oY/2 

Oo) - -er^^^|i^ln[(a;o-|H)/r],p = 2. (3.11) 

In the left plot of figured] we present the charge density as a function of the parameter aj:)/{2TT) in 
the D = A model with a single compact dimension of the length L £, . Note that for an untwisted scalar 
field this parameter is the flux measured in units of the flux quantum. For the chemical potential and 
for the length of the compact dimensions we have taken the values corresponding to ;U = 0.5m and 
mL£) = 0.5. The numbers near the curves correspond to the values of T/m. 



4 Current density 

Now we turn to the expectation value of the current density. As it can be easily seen, the components 
of the current density along the uncompactified dimensions vanish: (jV) = for r = 1, . . . ,p. By 
making use of Eq. p.2p and the expression (|2.26p of the Hadamard function, for the current density 
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Figure 1: The expectation values of the charge (left plot) and current (right plot) densities as functions 
of the parameter auj^'K for the Z? = 4 model with a single compact dimension and for /i = 0.5m, 
mLi) = 0.5. The numbers near the curves correspond to the values of T/m. 

along the r-th compact dimension we get: 

) = (2 W2+1 XI cosh(/xn/3) ^ rir sin(n,.d,.) 

E ^n^lVf +1K-1 Vn'rL'r+n^n, (4.1) 



X 

^9-1 



with r = p + 1, . . . , L) and, as before, the prime means that the term with n = should be taken 
with the weight 1/2. Note that, unlike to the case of the charge density, the current density does not 
vanish at zero temperature for a fixed value of the chemical potential. The zero temperature current 
density is given by the n = term in Eq. (|4.ip : 



The features of this current are discussed in detail in Appendix|Al For the model with a single compact 
dimension the general formula reduces to: 



(/) = ^ ''"^+1 E cosh(n^/3) ^ 
(2vr) 2 ^^^^ 



X sin(nrar)/ g+i (m^/re^L^ +n2/g2), (4.3) 

2 

An alternative expression of the current density is obtained by making use of the formula (j2.28p 
for the Hadamard function in Eq. (j3.2p : 



4gy-i ^ 
if) = if)o+ , Y.^osh{nfs(3) 
(2vr) 2 Lrn=l 

X J^(27rn, + a,)a;P;V£^(n/3^nJ. (4.4) 
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Prom Eqs. (|4.1|) and (|4.4|) it follows that, the current density along the r-th compact dimension is an 
odd periodic function of and an even periodic function oi ai, I ^ r, with the period equal to 27r. 
The current density is an even function of the chemical potential and it does not vanish in the limit 
of zero chemical potential. In the absence of uncompactified dimensions one has p = and from Eq. 
(j44]l we get 

if) = {f)o + TT-r Yl cosh(n;u/3) V (2^n, + a,) , (4.5) 

^^^^ n=l n, '^"^ 

where we have used f-i/2{x) = y/TT/2e~^. Here we assume that ojq > 0. In the case ujq = there is a 
zero mode and the contribution of this mode should be considered separately. 

In a way similar to that for the case of the charge density, we can see that in the limit when 
the length of the l-th compact dimension is large {I ^ r), the leading term obtained from Eq. ()4.4p 
coincides with the current density in the space with topology RP~^^ x {S^)'^~^ with the lengths of the 
compact dimensions Lp+i,. . . ,Lj_i,L;_|_i,. . . ,L£). For small values of Li, I ^ r, the behavior of the 
current density crucially depends whether di is zero or not. For di = the dominant contribution 
comes from the term with n; = and from the expression given above we can see that, to the leading 
order, Li (j^) coincides with the corresponding quantity in (D — l)-dimensional space with topology 
RP X {S^y^^ and with the lengths of the compact dimensions -Lp+i,. . . ,L;_i,Li+i,. . . ,L£). For a; 7^ 
and for small values of Li, the current density (j'") is exponentially suppressed. 

If ^ /3, the dominant contribution to the series over n in Eq. (14. 3p comes from large values 
of n ~ Lr//3. In this case we can replace the summation by the integration and the corresponding 
integral is evaluated by using the formula from Ref. \2Q] (assuming that < a;n;j_^). To the leading 
order we get 

) ^ (2n)(P+^)/^ ^ sm(n,a,) (a;„._^ - ^ ) 2 

^ ' 1=1 ni: 



xf{p+i)/2{nrLr^u;lr_^ - fj.^). (4.6) 

For a fixed value of Lr this formula gives the leading term in the high-temperature asymptotic for the 
current density. If in addition ^ Li, I ^ r, the dominant contribution comes from the term with 
Ur = 1, = 0, and the current density (j*") is suppressed by the factor e~ ""v '^or-A' ^ where 



In order to see the asymptotic behavior of the current density at low temperatures it is more 
convenient to use Eq. (|4.4p . Assuming that /3 (c^iq — |//|) ^ 1, the dominant contribution to the 
temperature dependent part comes from the mode with the smallest energy corresponding to = 
and one has 

P/2-1 -/3<^o + l/^l/3 



In this case the temperature corrections are exponentially small. 

For p ^ 2 the current density, defined by Eq. (j4.4p . is divergent in the limit |/i| — )• ujq. The 
corresponding asymptotic is found in a way similar to that for the case of the charge density. To the 
leading order we have 

(/)-^f^^(jo), (4.9) 

where the asymptotic expressions for (jo) for separate values of p are given in Eq. (13. lip . 

In the right plot of figure [T] we displayed the current density along the compact dimension as a 
function of d!£)/(27r) for the D = 4 model with a single compact dimension of the length corresponding 
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to mLi) = 0.5. The numbers near the curves are the values of T/m and for the chemical potential we 
have taken the value /i = 0.5m. 



5 Zeta function approach 

The expectation values of the charge and current densities can be evaluated directly from Eq. (13. ip 
by using zeta function techniques (see, for instance, Ref. |24]). First we consider the current density. 

5.1 Current density 

Substituting the expansion (|2.16p for the field operator and by making use the expression (j2.7p for the 
mode functions, for the current density along compact dimensions one finds the following expression 



p/3(wk-s/i) _ 1 



(5.1) 



with kr = (27rnr + cir) /L^ and r = p + 1, . . . ,D. The first term in the square brackets corresponds 
to the current density at zero temperature. The s = +/— terms are contribution coming from the 
particles/antiparticles. For the further transformations it is convenient to write Eq. (j5.ip in the form 



oo 

2e kr 



(f) = yl^y' n(Su^ cosh(n/3/i). (5.2) 



In the special case p = this formula is reduced to Eq. (j4.5p . In the representation (j5.2p . the zero 
temperature part corresponds to the n = term. The divergences are contained in this part only. 
The components of the current density along uncompact dimensions vanish. 
As the next step, in Eq. (|5.2p we use the integral representation 

—nf5uj n poo 

^ ' dse-^"'"-^"^"l^'\ (5.3) 







This allows us to write the expectation value of the current density in the form 



^ 2vr^ y ^ ^n/3,-n^/5V4.^ (5.4) 



Now we apply to the series over n the Poisson summation formula 

hoc 4-00 



^ g{na) = - 5(2vrn/«)> (5-5) 



a 

n=—oo 



where g{y) = j'^^dxe ^y^g{x). For the function corresponding to the series in Eq. ()5.4p one has 
g{y) = . After the integration over s we get the expression 

The current density defined by Eq. (15. 6p can be written as 

2 °° 

0") = 72 E (2^n, + a,)C.(s)U=i, (5.7) 

r 

' nr = — OO 
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with the partial zeta function 



Us) 



( 27rnD+i , . \ , 2 
I ^ h «/i 1 + m 



(5i 



where n!] 



/3y, y (27r)P 

(np+i, . . . , n^-i, nr-+i, . . . , n/j+i) and is given by Eq. (j2.1ip . Hence, in order to find 
the renormahzed value for the current density we need to have the analytic continuation of the zeta 
function (j5.8p at the point s = 1. 

The analytic continuation can be done in a way similar to that we have used in Ref. [9j for the 
zero temperature fermionic current. We first integrate over the momentum along the uncompactified 
dimensions: 



Us) 



T{s-p/2)Lr 
(47r)p/2r(s)y,/3 ^ 



K + 



( InriD 



+1 



/3 



+ z/x + m 



(5.9) 



Next, the direct application of the generalized Chowla-Selberg formula ^25j to the series in Eq. ()5.9 
leads to the following expression 



Us) 



T{s - D/2) 



(47r)«/2 T{s) (27r)^/2r(s) 
X 2^ cos(n^ • ag)/B_^(m^5n-(Lg)), 



(5.10) 



where L*^ 



and 



{Lp+l, . . . , Lr-l,Lr+l, . . . Ld+i), OLq — (Op+l, . . . , cir-l, Or+l, 

Ld+1 = OlD+l = 



QD+i), with 



(5.11) 



(5.12) 



The prime on the summation sign in Eq. (j5.10p means that the term = should be excluded from 
the sum and we use the notation 



1/2 



(5.13) 



for the vectors c = (ci, . . . , q) and b = (6i, . . . ,bi). Note that in Eq. (jS.lOp . cos(ng • cxg) can also be 
written as cosh(nD+i/x/3) H/ip+i.^r cos(niQ:i). 

The contribution of the second term on the right-hand side of Eq. (j5.10p to the current density 
is finite at the physical point. The analytic continuation is required for the part with the first term 
only. That is done, by applying the summation formula ()2.22p to the series over n^. The further 
transformations are similar to that we have used in deriving Eq. (I2.26P and we get 



T{s-D/2) 
(47r)^/2r(s) 



E 



22-s^D+3-2s^2 °o 

^^-y^^^^^y^^^nsin(na,0/i^_,(nL.m). 



(5.14) 



The right-hand side of Eq. (j5.14p is finite at the point s = 1. Now, substituting Eq. (j5.10p into Eq. 
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(|5.7p and using Eq. (|5.14|) . we find the following expression for the current density 

{f) = p 1^ ^ > n sm(nar)f D+i (nLrm) 

{27T)— ti 



(2vr)-L2 

cosing ■ dig) fD_^{mrgn'^^{^)). (5.15) 

Note that in the limit T — )• and Li ^ oo, I ^ r, the second term in the right-hand side of this 
formula vanishes. The first term presents the current density at zero temperature in the model with 
a single compact dimension (see Eq. (j4.2p for a special case p = D — 1). 

An alternative representation for the expectation value of the current density is obtained if we 
apply the formula (j3.4p to the series over rir in Eq. (|5.15p . Under the condition ^ m, this leads to 
the following expression 

= (27r)Wi)/2 z2 sin(n^a^) cosh(nD+i///3) 



X cos{nl_^ ■ al_-^)fD±i{m^ gl^(Lg) + n%^^f3'^), (5.16) 

where a^^i = {dp+i, . . . ,dr-i,dr+i, ■ ■ ■ do), Lg = {Lp+i, . . . , Ld), and gn^(Lg) is defined by Eq. 
(|5.13p . In particular, for a massless field and for zero chemical potential, fJ- = 0, from (|5.16p we get 



r((ij + i)/2) ^ ^ 

(j ) = 2eL,. — ^^^^^^^ — 2^ nr.sm(nrar 



oo 



n,.=l 

V cos(n^_i • ag_i) 

(L,) + n2/32](i5+i)/2- ^ • ) 



X 

1-5 



The equivalence of two representations for the current density, Eqs. (|4.ip and (j5.16p . can be seen 
by using the relation 



J2 cos(n • a)/.(cJ&2 + ^' «2^2) ^ ^ «;r"'/.-V2(fe^n), (5.18) 

n n 

where n = (ni, . . . ,ni), a = {ai, . . . , a;), and w'^ = ^[=i(27rnj + ai)'^/af + c^. This relation has been 
proved in Ref. |17] by using the Poisson's resummation formula. Note that the formula (13. 4p is a 
special case of Eq. (|5.18p . 

An expression for the current density, convenient for the discussion of the high-temperature limit, is 
obtained from Eq. (j5.16p . by applying to the series over n/j+i the formula (|3.4p under the assumption 

^ m. This leads to the following expression 

A L 

^j'''> = (n \D/2R sin("'rar) Y cos(n^_i • a^-i) 

V^^'' P nr=l n- 



[(27rnD+i//3 + ifif + mY^^fD/2{9n, (L,) V(27rnD+i//? + ifif + m^). (5.19) 
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At high temperatures the dominant contribution comes from ud^i = term and to the leading order 
we have 

^ 7?W2 £ nrsm{nrdr) cos(n^_i • dcq^i) 

X (m2 -;x2)^/Vd/2K(L,)x^^^^). (5.20) 

The corrections to this leading term are exponentially small. The equivalence of two representations, 
Eqs. (j4.6p and (j5.20p . for the leading order term can be seen by using the relation (j5.18p . 

5.2 Charge density 

Now we turn to the evaluation of the charge density by using the zeta function approach. Similar to 
the case of Eq. ()5.ip . we have the following mode sum 

= WFV, ^ £ e/3(-.-M)-l- (^-^^^ 

The zero temperature part in the charge density vanishes due to the cancellation between the con- 
tributions from the virtual particles and antiparticles. The corresponding contributions to the finite 
temperature part have opposite signs due to the opposite signs of the charge for particles and antipar- 
ticles. Introducing the expectation values for the numbers of the particles and antiparticles (per unit 
volume of the uncompactified subspace), 

(^^) = 7^E ,,(J.)_1 ^ (5-22) 



k 



the charge density is written as (jo) = e {N^ — N^) /Vg. In Eq. (j5.22p . the upper/lower sign corre- 
sponds to particles/antiparticles. Note that in the current density the contributions from particles 
and antiparticles have the same sign (see Eq. (|5.ip ). This is due to the fact that, though the charges 
have opposite signs, the opposite signs have the velocities as well, fr"*"^ = k-r/u; for particles and 
vi ^ = —kr/tiJ for antiparticles (see the phases in the expression (j2.7p for the mode functions). The 
expression for {N±) is obtained from Eq. ()3.5p by the replacement 2e sinh(n/i/3)/ — t- e^"'^^. 
The expression (j5.2ip for the charge density may be written in the form 

oo 

= (2^ E E ^"'^"^^ ^-h(n/3/.). (5.23) 

^ ' ''^ k n=l 

For the further transformation of this expression we use the relation 

^npuj ^ Jq ) cosh(n/3^). (5.24) 

As a result, the expectation value of the charge density is presented in the form 

(jo> = [fi- dfidpfij ^^^__cosh(n/3^). (5.25) 

Substituting Eq. ()5.3p . by the transformations similar to that we have used in the case of the 
current density, one finds 



Oo) 



2e (^fi- £ dfidp(3^ C{i 



(5.26) 
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where the corresponding zeta function is defined as 



Cis) 



1 



dk. 



27rn 



(5.27) 



with cjk defined by Eq. (12. lip . After the integration over over the momentum along uncompact 
dimensions, the function ()5.27p is written in the form 



as) = ^(izjVHL V 



i=p+i \ J 



+ m 



(5.28) 



where rig+i = (n^+i, . . . ^nD+i) and L^^i, ao+i are defined by Eq. (j5.1ip . The application of the 
generalized Chowla-Selberg formula [25] to Eq. (j5.28p gives 



C(^) 



^^^+i-2. r(g-p + i_)/2) ^ 21 



m 



D+1-2S 



^9+1 



(47r)(^+i)/2r(s) (27r)(-D+i)/2r(s) 



X > COS(nq+i • Q;q+i)/D±l_^ (^my'c/^^(Lg) +n^_^^ 



(5.29) 



with Lg+i = (-Lp+i, . . . ,L£)_|_i) and Qg+i = ((ip+i, . . . ,aD+i)- The prime on the summation sign in 
Eq. (|5.29p means that the term with ni =0, / = p + 1, . . . , -D + 1, should be excluded from the sum. 
Substituting Eq. (j5.29p into Eq. (j5.26p . for the charge density one finds the expression 



Oo) 



re sinh(^/3n) 



4erei^+i/3 
(27r)(^+i)/2 



^cos(ng • aq)fD±i {m^ gl^{l^q) +rfip'^^ . 



n=l 



(5.30) 



Note that the first term in the right-hand side of Eq. (I5.29P does not depend on temperature and the 
corresponding contribution in Eq. ()5.25p vanishes. This expression for the charge density is valid for 
the region ^ m. The equivalence of the representations (13. 5p and (j5.30p in this region is proved 
by using the formula (j5.18p . In Eq. (j5.30p . the term with re; = 0, / = p + 1, . . . , D, coincides with the 
corresponding charge density in Minkowski spacetime {p = D, q = 0) given by Eq. (13. 6p . Note that 
by the replacement 2e sinh(re/i/3)/V^ — ?• e^"'^^ in Eq. ()5.30p . we can obtain the corresponding formula 
for {N±). 

An alternative expression for the charge density, convenient for the investigation of the high- 
temperature limit, is obtained from Eq. (j5.30p if we first separate the part corresponding to (jo)(m) 
and then apply to the series over n in the remained part the formula (13. 4p . This leads to the following 
expression: 



Oo) 



Oo) 



2ieT 



(M) 



:[(27rrer + i^f + m"^]^-^ fD_-^{gn,{Lg)^y {2TTnT + ifi)^ + m^). 



(5.31) 



As before, the prime means that the term with rii = 0, I = p -\- 1, . . . , D, should be excluded from 
the sum. At high temperatures the dominant contribution to the second term in the right-hand side 
comes from the term with re = 0: 



Oo)m + 



(2^)^/2 



E' 



X (m — ^J■ ) 



2^f 



-1 



/D_l(5n,(Lg)\/"l2 - /x2) 



(5.32) 



15 



The higher order corrections to this asymptotic expression are exponentially small. Hence, similar to 
the case of the current density, the topological part of the charge density is a linear function of the 
temperature in the high-temperature limit. 

In order to find the asymptotic expression for the part (jo)(m) high temperatures, we use the 
integral representation 

Substituting this into Eq. p.6p and changing the orders of integration and summation, the summation 
is done explicitly and to the leading order we get 

Oo>(M) « 2ef^T^-'^^^^^^j±^Cn{D - 1), (5.34) 

where Cr(3^) is the Riemann zeta function. This result has been obtained in Ref. |23j . As it is seen 
from Eq. (j5.32p . for D > 2 the Minkowskian part dominates in the high temperature limit and one 

has {jo) f« (io)(M)- 

6 Bose-Einstein condensation 

In this section we consider the application of the formulas given before for the investigation of the 
Bose-Einstein condensation (BEC). This phenomena for a relativistic Bose gas of scalar particles in 
topologically trivial flat spacetime has been discussed in Refs. [22l [231 [26] . The investigation of the 
critical behavior of an ideal Bose gas confined to the background geometry of a static Einstein universe 
is given in Refs. |27j for scalar and vector fields. BEC in higher dimensional spacetime with as 
a compact subspace has been considered in Ref. [28]. The case of an ultrastatic (3+l)-dimensional 
manifold with a hyperbolic spatial part is analyzed in Ref. [29]. The background geometry of closed 
Robertson- Walker spacetime is discussed in Ref. [30]. Thermodynamics of ideal boson and fermion 
gases in anti-de Sitter spacetime and in the static Taub universe have been considered in Refs. [3T|[32]. 
In the high-temperature limit, BEC in a general background has been discussed in Refs. \33 \ 134 1 [35]. 
Recently, BEC on product manifolds, when the gas of bosons is confined by anisotropic harmonic 
oscillator potential, has been investigated in Ref. [36J. 

Note that in the literature two criteria have been considered for BEC (see, for instance, the 
discussion in Refs. [37J). In the first one the existence of critical temperature Tc > is assumed for 
which the chemical potential becomes equal to the single particle ground state energy. The derivative 
drfJ' for a fixed value of a conserved charge is discontinuous at the critical temperature and the 
condensation corresponds to a phase transition. By the second criterion, one assumes the existence 
of a finite fraction of particle density in the ground state and in states in its neighborhood at T > 0. 
In this case the presence of a phase transition is not required and thermodynamical functions can be 
continuous. In particular, in Ref. [23j it has been shown that for massive particles there is no BEC in 
dimensions D ^ 2 if one follows the first criterion. 

In the discussion of previous sections we have considered the charge and current densities as 
functions of the temperature, chemical potential and the lengths of compact dimensions. From the 
physical point of view it is more important to consider the behavior of the system for a fixed value 
of the charge. We will denote by Q the charge per unit volume of the uncompactified subspace, 
Q = Vq {jo)- From Eq. p.Sp for this quantity one has 

4/3 ^ 

For a fixed value of the charge, this relation implicitly determines the chemical potential as a function 
of the temperature, lengths of the compact dimensions and the charge. 



16 



For high temperatures the chemical potential determined from Eq. (|6.1|) tends to zero. Hence, 
at high temperatures we always have solution with |/i| < ojq. The further behavior of the function 
//(T) with decreasing temperature is essentially different in the cases p > 2 and p ^ 2. For p > 2 the 
expression (|6.ip is finite in the limit — >• ujq. We denote by Tc the temperature at which one has 
|/i(rc)| = ujQ for a fixed value of the charge. It is the critical temperature for BEC. The formula 



IQI 



4|e|/3e 
(2vr)^ 



- n sinh(nLi;o/3c) 



n=l 



fp+l 
2 



(n/3cWn 



(6.2) 



with /3c = 1/Tc, determines the critical temperature as a function of the charge, of the lengths of the 
compact dimensions, and of the parameters a;. 

Simple asymptotic formulas for the critical temperature are obtained for low and high temper- 
atures. At low temperatures, ujqPc ^ Ij the dominant contribution in Eq. (j6.2p comes from the 
mode with the smallest energy corresponding to ni = 0. By using the asymptotic expression for the 
Macdonald function for large values of the argument, from Eq. (j6.2p to the leading order one finds 



2tt 



\Q/e\ 

Cr(p/2) 



2/p 



(6.3) 



This regime is realized for values of the charge corresponding to \Q/e\ <C ujq. At high temperatures. 



by taking into account that to the leading order (jo) 
(|5.34p . one finds 

r(^+i)/2|Q/e|K- 



Oo)(M) using the asymptotic expression 



2wor((Z) + l)/2)CRiD - 1) 



1/{D-1) 



(6.4) 



This asymptotic corresponds to \Q/e\ 3> VqUj^j'. In figure [2] we display the critical temperature as a 
function of the charge density in the D = 4 model with a single compact dimension (p = 3, g = 1). The 
graphs are plotted for mLo = 0.5 and for different values of aD/i2iT) (numbers near the curves). As it 
is seen, for fixed lengths of the compact dimensions, the critical temperature for the phase transition 
can be controlled by tuning the value for the gauge potential. 




Figure 2: The critical temperature as a function of the charge density in the D = 4 model with a 
single compact dimension. The graphs are plotted for Lom = 0.5 and for different values of and 
the numbers near the curves correspond to the values of the parameter a£)/{2T:). 

At temperatures T < Tc, the equation (j6.ip has no solutions with < ojq. The consideration in 
this region of temperature is similar to the standard one for the BEC in topologically trivial spaces. 
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We note that the expression (|6.ip does not include the charge corresponding to the states with kp = 0. 
At temperatures T < the expression (|6.ip with |//| = determines the charge corresponding to 
the states with kp / 0. We denote this charge by Qi: 



4e/3sgn(^) m p+i / t a 

(27r) ^ 



3.5) 



n=l 



For the charge corresponding to the Bose-Einstein condensate at the state kp = one has Qc = Q — Qi. 
This charge vanishes at T = T^- At low temperatures, by making use of Eq. (|6.3p . for the corresponding 
charges below the critical temperature one finds 



Qi = Q{T/T,y'\ Q, = q\i- {T/TcY'^' 

For high temperatures we use the asymptotic formula (|6.4p with the results: 

Qi = Q{T/T,)''-\ Q, = Q[l- {T/T,f-'] . 



(6.6) 



(6.7) 

In particular, Eq. (j6.7p coincides with the corresponding result in Minkowski spacetime [23] . In figure 
[3l we plot the chemical potential as a function of the temperature in the D = 4 model with a single 
compact dimension {p = S, q = 1). The left and right plots correspond to od = and an = 7r/2 
respectively. For the length of the compact dimension we have taken the value corresponding to 
mLf) = 0.5 and the numbers near the curves correspond to the values of the parameter m}~^Q/e. 





Figure 3: The chemical potential as a function of the temperature in the L> = 4 model with a single 
compact dimension for = 0.5. The left and right plots correspond to old = and a/j = 7r/2 

respectively. The numbers near the curves correspond to the values of the parameter m}-~^Q/e. 

Note that for T < the scalar field acquires a nonzero ground-state expectation value ipc- The 
latter can be found in a way similar to that used in Ref. [33] (see also Refs. |22] 123] ): 

with Q and Qi given by Eqs. (16. ip and ()6.5p . 

Having the chemical potential from Eq. (j6.ip for T > Tc and taking = ujq for T < Tc, we 
can evaluate the current density for a fixed value of the charge by using Eq. (14. 4p . Note that at 
high temperatures the chemical potential tends to zero and the leading term in the corresponding 
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asymptotic expansion for the current density is obtained from Eq. (|5.20p with ^ = 0. For T < T^, the 
formula (j4.4p gives a part of the current density due to the excited states only. In addition to this, 
there is a contribution due to the condensate, given by the expression 



^rQc 
LrUJ^Vq 



(6.9) 



with r = p + 1, . . . , D and la^l ^ vr. For the expectation value of the total current density one has 

{f)t=fc + {f)- 

In figure[3]we plot the expectation values of the total current density (left plot, full curve) and of the 
particle and antiparticle numbers in the states with kp ^ 0, (A'^-i-) (right plot), versus the temperature 
for mLi) = 0.5, ud = ''r/2 and Q/e = 0.5m^~^ in the D = A model with a single compact dimension. 
For the corresponding critical temperature from Eq. (16.20 one finds ~ 0.63m. In the left plot we 
have separately presented the contributions to the current density coming from particles (dot-dashed 
line) and antiparticles (large dashed line). The linear dependence at high temperatures is clearly seen. 
On the left panel we have also plotted the separate contributions to the current density from the 
excited states, {j^) /(em^) (dashed line), and from the condensate, j^/{em^) (dotted line). Note 
that for the current density given by Eq. ()4.2p we have (j^)q = 2.32em^. The total current and 
its first derivative with respect of the temperature are continuous functions at the point of the phase 
transition. The latter is not the case for the separate parts coming from the condensate and from the 
excited states. 





Figure 4: The expectation values of the current density along the compact dimension (left plot) and 
of the particle/antiparticle numbers (right plot) as functions of the temperature in the D = 4 model 
for mLr) = 0.5, an = 7'"/2 and Q/e = 0.5m^~^. 

For p ^ 2 the charge defined by Eq. (jG.ip diverges in the limit — )• and for finite value of 
charge density the point |//| = m cannot be reached. The corresponding asymptotic expressions for 
the charge and current densities are given by Eqs. p. lip and ()4.9p . Thus, for p ^ 2 there is no BEC 
by the first criterion given above. In particular, this is the case in the model with compact space 
corresponding to p = 0. This result for spaces of finite volume in general has been obtained in Ref. 
|35] . As before, for p ^ 2 and for a fixed value of the charge Q, the dependence of the chemical 
potential on the temperature is determined by Eq. ()6.ip . In the limit /3 — ?• oo (low temperatures) and 
for a fixed value of |/i| ^ loq, the expression on the right-hand side of Eq. (|6.ip tends to zero. From 
here we conclude that for a fixed value of Q we should have \fi\ — )• cjq for /3 — )• oo. The corresponding 
asymptotic behavior is found in a way similar to that we have used for Eq. (|3.10p and is given by 
the same expression. Solving with respect to the chemical potential, we find the following asymptotic 
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expressions 



l/^l 



Wo 



1^0 \ ~-p 
271, 



Wo — T exp 



\Q\ \ 2) 

IQI f2^V"' 
|e| \ujqT J 



2 

2-p 



p = 0,l, 



p = 2, 



(6.10) 



in the limit T — )• 0. In figure we have plotted the chemical potential versus temperature in the 
D = 3 model with a single compact dimension [p = 2, q = 1) for ajj = and for a fixed value of the 
charge density corresponding to m}^^ (jq) /e = 0.5. The numbers near the curves correspond to the 
values of the parameter Ldtii. As it is seen, though for a finite value oi Ln the derivative dTfJ-iT) 
is a continuous function, in the limit L^) — )• oo it tends to the corresponding function in the D = 3 
model with trivial topology {p = 3, q = 0) for which the function dTfJ-iT) is discontinuous at the point 

T = Tr. 




Figure 5: The chemical potential as a function of the temperature in the D = 3 model with a single 
compact dimension for a/? = and m}^^ (jo) /e = 0.5. The numbers near the curves correspond to 
the values of the parameter Lom. 



7 Conclusion 

In the present paper we have investigated the finite temperature expectation values of the charge and 
current densities for a complex scalar field, induced by nontrivial spatial topology. As an example for 
the latter we have considered a flat spacetime with an arbitrary number of toroidally compactified 
dimensions. This allowed us to escape the problems related to the curvature and to extract pure 
topological effects. The periodicity conditions along compact dimensions are taken in the form (j2.2p 
with general constant phases. As special cases the latter includes the periodicity conditions for 
untwisted and twisted fields. In addition, we have assumed the presence of a constant gauge field. By 
performing a gauge transformation, the gauge field is excluded from the field equation. However, this 
leads to the shift in the phases appearing in the periodicity conditions given by Eq. (12. 5p . 

In the evaluation of the expectation values for the charge and current densities we have used 
two different approaches which allowed us to obtain alternative representations for the corresponding 
expectation values. In the first approach we evaluated the thermal Hadamard function by using the 
Abel-Plana-type summation formula for the series over the momentum along a compact dimension. 
The corresponding expression is given by Eq. ()2.26p . The = term in that formula corresponds to 
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the Hadamard function for the topology R^^^ x (5^)"^^^ and, hence, the n,. 7^ part is the change in 
the Hadamard function due to the compactification of the r-th direction. An alternative representation 
for the Hadamard function is given by Eq. (I2.28|) . 

Given the Hadamard function, the expectation values of the charge and current densities are 
evaluated by making use of Eq. (j3.2p . The charge density is given by two equivalent representations, 
Eqs. (j3.3p and (j3.5p . The explicit information contained in Eq. (|3.3p is more detailed. The term 
with nr = in this representation corresponds to the charge density for the topology i?^"^^ x (S^)'^~^ 
with the lengths of the compact dimensions (Lp^i, . . . , Lr-i, L^+i, . . . ,Ld) and the contribution of 
the terms with rir 7^ is the change in the current density induced by the compactification of the r-th 
dimension. The charge density is an even periodic function of the phases ai with the period equal to 
27r. The sign of the ratio (jo) /e coincides with the sign of the chemical potential. If the length of the 
Z-th compact dimension is small compared with the other length scales and Li <C /?, the behavior of 
the charge density is essentially different with dependence whether the parameter ai is zero or not. For 
ai = 0, to the leading order, Li (jo) coincides with the charge density in (D — l)-dimensional space with 
topology RP X (S^y~^ and with the lengths of the compact dimensions Lp+i,. . . ,Li_i,Li^i,. . . ,L£). For 
ai the charge density is suppressed by the factor e"'"''^/^'. At low temperatures and for a fixed 
value of |//| < cuq, the charge density is suppressed by the factor e^^'^"^!'^!)/-^. For a fixed temperature 
and in the limit — t- ujq the charge density is finite for p > 2 and it diverges for p ^ 2. The 
corresponding asymptotic behavior in the latter case is given by Eq. (j3.1ip . In the high-temperature 
limit and for D > 2 the Minkowskian part dominates in the charge density with the leading term 
given by Eq. (j5.34p . In the same limit, the topological part of the charge density is a linear function 
of the temperature. 

For the expectation value of the current density along the r-th compact dimension we have derived 
representations given by Eqs. (j4.ip and (j4.4p . The components along uncompactified dimensions 
vanish. The current density along the r-th compact dimension is an odd periodic function of Or 
and an even periodic function of a/, / 7^ r, with the period equal to 2tt. The current density is an 
even function of the chemical potential. Unlike to the case of the charge density, the current density 
does not vanish at zero temperature for a fixed value of the chemical potential. The corresponding 
expression is given by Eq. (14. 2p and the properties are discussed in Appendix |Al For small values of 
Li, I ^ r, and for ai = 0, to the leading order, the quantity Li {j^') coincides with the r-th component 
of the current density in {D — l)-dimensional space with topology x and with the lengths 

of the compact dimensions Lp+i,. . . ,L;_i,L;_|_i,. . . ,L£). For a; 7^ and for small values of Li, I 7^ r, 
the current density (j*") is exponentially suppressed. At a fixed temperature and for p ^2 the current 
density is divergent in the limit — )• coq. The leading term in the corresponding asymptotic expansion 
is related to the charge density by Eq. (j4.9p . For a fixed value of the chemical potential |//| < coq and at 
low temperatures the finite temperature corrections are given by Eq. ()4.8p and they are exponentially 
small. In the limit of high temperatures, the current density is a linear function of the temperature. 

In Section m we have derived alternative representations for the expectation values of the charge 
and current densities by using the zeta function approach. In both cases, by applying to the corre- 
sponding zeta functions the generalized Chowla-Selberg formula, Eqs. (15.16p . (15.19P and Eqs. (15.30p . 
()5.3ip are obtained for the current and charge densities respectively. At high temperatures, the lead- 
ing term in the asymptotic expansion of the current density is given by Eq. (I5.20p with the linear 
dependence on the temperature and the next corrections are exponentially small. For the charge den- 
sity, for D > 2 the leading term in the high-temperature expansion coincides with the corresponding 
charge density in (D + l)-dimensional Minkowskian spacetime. The leading term in the correction 
induced by nontrivial topology linearly depends on the temperature and the following corrections are 
exponentially suppressed. 

The Bose-Einstein condensation is discussed in section [6j For a fixed value of the charge, the 
relation ()6.ip determines the chemical potential as a function of the temperature, of the lengths of 
compact directions and of the phases in the periodicity conditions. For high temperatures the chemical 
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potential tends to zero. With decreasing temperature the chemical potential increases and for p > 2 
one has = at some finite temperature T = Tc. The critical temperature for BEC, Tc, is 

determined by Eq. (16. 2p . Simple expressions are obtained for low and high temperatures, Eqs. ()6.3p 
and ()6.4p . respectively. At temperatures T <Tc one has |^| = a;o and Eq. (j6.ip determines the charge 
corresponding to the states with kp 7^ and the remained charge corresponds to the charge of the 
condensate. At low and high temperatures the charges are given by simple expressions (|6.6p and (j6.7p . 
Similar to the charge density, for T < Tc the current density is the sum of two parts. The first one is 
the contribution of excited states and is given by Eq. (j4.4p with |//| = loq. The second part is due to 
the condensate and it is presented by Eq. (16. 9p . The total current and its first derivative with respect 
to the temperature are continuous functions at the critical temperature. For p ^ 2, the point = ujq 
cannot be reached for finite value of charge density. For a fixed value of the charge, we have — )• 
in the limit T — >• 0. The corresponding asymptotic behavior is given by Eq. (I6.10p . In this case the 
thermodynamical functions are continuous and there is no phase transition at finite temperature. 
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A Vacuum expectation value of the current density 

In this Appendix we give some properties of the zero temperature current density given by Eq. (14. 2p . 
An alternative expression is obtained from Eq. (I5.16|) taking the limit /? — )> 00. In this limit the term 
with n = survives only and we get: 

)o = (2^)(D+i)/2 2^^^rsm(n,a,) 

^ cos(n^_i • OLq-\)f D±i{jng^^{J^q)). (A.l) 



X 

Let us consider the behavior of the zero temperature current density in some limiting cases. First we 
consider the limit when the length of the r-th compact dimension, L^, is much larger than the other 
length scales. The behavior of the current density in this limit crucially depends whether defined 
by (14. 7p . is zero or not. In the first case, which is realized for = 0, Z 7^ r, and m = 0, the dominant 
contribution in Eq. (j4.2p for large values of Ly. comes from the modes with ni = 0, I ^ r, for which 
Wn^ i = ojQr = 0. The corresponding expression is obtained from Eq. (14. 2p taking the limit Wn^_^ 
and to the leading order we have 

_ 2eT{p/2 + 1) ^ sin(n^a^) 

For ujqj. 7^ and for large values of L.^, the main contribution to the zero temperature current density 
comes from the mode = 1, n/ = 0, / 7^ r, and from Eq. (j4.2p one finds 

2eK-isin(a,)4^+^)/' , 

In this case we have an exponential suppression. 

Now we discuss the asymptotic of the current density for small values of Lj-. In this limit it is 
more convenient to use Eq. (jA.ip . First we separate the term n; = 0, Z 7^ r, in Eq. (jA.ip and use the 
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asymptotic expression of the Macdonald function for small values of the argument. For the remained 
part in Eq. (jA.ip . the dominant contribution comes from large values of rir and, to the leading 
order, we replace the summation over rir by the integration. The corresponding integral involving the 
Macdonald function is evaluated by using the formula from Ref. |20j. In this way it can be seen that 
the contribution of the mode with a given ng_^ is suppressed by the factor exp{— gr ^ a'^ / L"^ + m?), 

where gr = 'Yld=p+i "-F-^r ^ result, we see that the dominant contribution to the current density 
is due to the modes with ni = 0, I ^ r, and to the leading order we get 

. ^ 2er((Zj + l)/2) ^ sin(n,a,) 

' n,. = 1 ' 

This leading term does not depend on the mass and on the lengths of the other compact dimensions. 
As it is seen from Eq. (jA.ip . the expression in the right-hand side of Eq. ()A.4p coincides with the 
current density for a massless scalar field in the space with topology R^^^ x S^. 
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